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Correlation functions (CFs) of the classical Heisenberg antiferromagnet on the pyrochlore 
lattice are studied by solving exactly the infinite-component spin-vector model. As in many 
Fully Frustrated Lattices, the constraint due to the minimization of the energy and the particular 
structure based on corner sharing tetrahedra both contribute to the creation of local degrees 
of freedom. The resulting degeneracy destroys any magnetic order at all temperature and we 
obtain no sign of criticality, even at T = 0. Calculated neutron scattering cross sections have 
their maxima beyond the first Brillouin Zone and reproduce experimental results obtained on 
Y(Sc)Mn2 and CsCrNiFg as well as theoretical predictions previously obtained by classical 
Monte Carlo simulations. Evidences for thermal and spatial decoupling of the magnetic modes 
are found so that the magnetic fluctuations in this system can be approximated by S{q, T) « 
f{l)h{T). 

PACS numbers: 75.10.Hk, 75.50.Ee, 75.40.Cx, 75.40.-s 



I. INTRODUCTION 

Motivated by the search and understanding of new 
quantum or classical disordered ground states in mag- 
netic systems, recent years have seen a renewal of interest 
in the properties of frustrated Heisenberg antiferromag- 
nets. Special attention focused on intrinsically frustrated 
models, where frustration is mainly coming from the con- 
nectivity of the underlying lattice and not only from com- 
peting interactions. From this point of view, two lattices 
distinguished themselves : the kagome lattice (made of 
corner sharing triangles) and the pyrochlore lattice (made 
of corner sharing tetrahedra) . Experimental and the- 
oretical studies related to both lattices have shown 
that materials as well as models exhibit unconvention- 
nal magnetic properties, involving noncollinear or incom- 
mensurate orderings, apparently broken ergodicity with- 
out structural disorder and in general, deviations from 
canonical behaviors. For most geometrically frustrated 
systems, the absence of ordering at finite temperature or 
the tendancy to a reduction of the real interactions seems 
to be a common characteristic. An experimental signa- 
ture of this property is a persistent paramagnetic state 
down to low temperature without any sign of ordering B] 
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and sometimes a transition to a frozen state at low tem- 
perature {Tp « |Ocvk|), with an unconventional 
slowing down of the dynamics H, ^. The theoreti- 
cal understanding of these properties is in progress but 
it seems that several properties, expected to be generic, 
are not systematically reproduced and depend on the ma- 
terial or the model. One point is at least widely shared: 
due to frustration, short range correlations develop at a 
temperature much lower than the interactions. This phe- 
nomenon has been proven by mean field analysis Q and 
classical monte-carlo simulations ||, |l^ . On the other 
hand, there are some properties much more dependent 
on the studied system. A major theme of study has been 
the effect of "order by disorder" previously quoted by 
Villain and co-workers and also by Shender It 
was noticed that switching on any source of fluctuation 
(like thermal or quantum fluctuations) has a tendancy to 
lift the degeneracy between the low lying classical ground 
states. Indeed, thermal fluctuations give rise to entropic 
selection in the kagome case Q but this kind of selec- 
tion was shown to be fluctuation-like dependent , and 
sometimes absent [16, |l^. Other types of perturbation 
have been taken into account, such as anisotropy [[l8[ , 
TqI , longer range interaction s |]7| , ]2C|] 
Til El |2^, |24|, E5|, ^ 



structural disorder 
and quantum fluctuations 

leading in general to a magnetic mode selection. 

The pyrochlore lattice is a three dimensional arrange- 
ment of corner sharing tetrahedra as depicted in Fig. ^. 
According to Lacorre definition of frustration it is 
probably the most frustrated structure. Experimentally, 
several families of compounds are known to crystallize 
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with this geometry. The first one is an oxide family with 
the general formula A^+B^+Ot (Refs. |, |9[ |ol) (A is 
a Rare Earth ion and B a transition metal). The sec- 
ond one concerns spinels of general formula AB20i or 
A2B2O4 (Refs. |l[ H ll) (A is a Rare Earth ion and 
B a transition metal). The third one concerns fluorides 
of general formula AB^+C^+Fg (Refs. |§ ||) (A is an 
alkali metal, B and C are transition metals). The last 
one is the intermetallic Laves Phases compounds of gen- 
eral formula AB2 (Refs. |3^, (A is a Rare Earth ion 
and B is a transition metal). Strictly speaking, the only 
candidate that can be rigorously described by a localized, 
uniform lattice model is the first family, as this system is 
well ordered with a uniform distribution of magnetically 
localized ions. The second and third sets often possess 
positional or chemical disorder |31, Q that should be 
taken into account in the model and the last one is bet- 
ter described within an itinerant model |38[ . 

In this paper, we will focus only on the isotropic classi- 
cal Heisenberg model on the pyrochlore lattice, neglecting 
all other possible contributions that might be necessary 
for an understanding of a peculiar material. This choice 
is supported by our hope to deduce intrinsic properties 
which come exclusively from geometrical structure. Our 
goal is to compute correlation functions using an approx- 
imation of the Heisenberg model, the infinite component 
spin- vector model. 



II. LATTICE STRUCTURE AND MODEL 

The "next simplest" approximation for classical spin 
systems, which follows the Mean Field Approximation 
(MFA), consists in generalizing the Heisenberg Hamilto- 
nian for the D-component spin vectors: 
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(1) 



and taking the limit D ^ 00. In this limit, it is possible 
to take into account exactly the effect of the Goldstone 
modes or would be Goldstone modes, which are the most 
important modes in frustrated magnets. We will only 
give the outline of the calculation here and give in a fur- 
ther article the details of the technique. Nevertheless, 
readers who are interested in the full derivation can find 
a description of the method in Ref. where the antifcr- 
romagnetic model on the kagome lattice was considered. 

The pyrochlore lattice shown in Fig. |^ consists of 
corner-sharing tetrahedra. Each node of the underly- 
ing Bravais lattice (face centered cubic) is numbered by 
i = 1,... ,N and corresponds to a tetrahedron. Each 
site of elementary tetrahedra is labeled by the index 
I = 1,2,3,4. 

It is convenient to put the Hamiltonian (^ into a diag- 
onal form. First, one goes to the Fourier representation 




FIG. 1: Structure of the pyrochlore lattice. The elemen- 
tary tetrahedra are labeled by i = 1, . . . , A'', and the sites on 
tetrahedra (the sublattices) are labeled hy I — 1, 2, 3, 4. 
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where the wave vector q belongs to the I''' Brillouin zone 
of the underlying face centered cubic lattice (see Fig. ||). 
The Fourier-transformed Hamiltonian reads 



n - — Vy"'s' -s'' 



(3) 



where the interaction matrix is given by 



Vq = 2J 
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.qx - g^^ 

e = cos( ) 



,qx+qy. „ qy-q^ 
c = cos( — -) / = cos(-S^— ) 



At the second stage, the Hamiltonian (^ is finally diag- 
onalized to the form 



(5) 



where — 2Jt'„(q) are the eigenvalues of the matrix 
taken with a negative sign, 

^3,4 = -1 ± v/l + Q- (6) 



1^1,2 = 1, 
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FIG. 2: Brillouin zone of the pyrochlore lattice. Letters F, 
X, W, U and L are high symmetry points using the usual 
crystallographic conventions. 



FIG. 3: Reduced eigenvalues of the interaction matrix, 
Vn{<l) = V'q"/(2J) of Eq. (|). They are plotted such that 
they cover the high symmetry axis of the face centered cubic 
Brillouin Zone (see Fig. H). 



where 



•3=2^ 2 2 ^ 

cos( ^)+cos( 

cos( — - — ) + cos(^- )) 



(7) 



2 ' '2 
The diagonaUzing transformation has the explicit form 



U-,\q)V^'Ul,rA^)^K^nn', 



(8) 



where the summation over the repeated indices is implied 
and U is the real unitary matrix, = , i.e., U^i^ = 

Uln- 

The Fourier components of the spins Sq and the col- 
lective spin variables cr" are related by 



Sq = C//n(q)crq, < = s^C/,„(q). 



(9) 



The largest dispersionless eigenvalues i'i^2 of the inter- 
action matrix [see Eq. (||) and Fig. |j are a manifestation 
of frustration in the system which precludes an extended 
short-range order even in the limit T 0. Since i>i^2 
are independent of q, 1/2 of all spin degrees of freedom 
are local and can rotate freely. The eigenvalue which 
becomes degenerate with z/i_2 in the limit g — > is re- 
lated, to the usual would be Goldstone mode destroying 
the long-range order in low-dimensional magnets with a 
continuous symmetry whereas the eigenvalue is posi- 
tioned, in the long-wavelength region, much lower than 
the first two ones, and it is tempting to call it the "opti- 
cal" eigenvalue. 



III. RESULTS 

Following the derivation introduced in Ref. ^ we 
compute all CF's at any temperature and in particular, 
the static magnetic neutron scattering cross section : 



da 



(10) 



where Q is a vector of the reciprocal space and r and 
r design spin belonging to the pyrochlore lattice, i.e., 
is a spin located on site i of the Bravais lattice 
in the position I of the unit cell. Two peculiar planes 
of the reciprocal space have been investigated as shown 
in Fig. ^ and Fig.^. The major feature is that there is 
almost no magnetic signal in the first Brillouin zone, all 
the weight being recovered in soft cusps outside. In both 
planes, we observe patterns very similar to the one ob- 
tained by Monte Carlo simulations in Ref. ^ or by neu- 
tron diffraction on Y{Sc)AIn2 (2j. We also note that the 
cross section in the ([hhO],[001]) plane is very close to the 
one previously calculated in the kagome case by the same 
method ]40[ |. The main implications of these results are 
that the quantum or classical nature of the spins in the 
pyrochlore lattice does not seem to be of relevance for un- 
derstanding the liquid nature in terms of spin-spin corre- 
lations. Indeed, different methods and measures, applied 
to differents spin species (classical or quantum) , give the 
same magnetic responses. Besides, the very strong re- 
semblance between the neutron diffraction pattern in the 
([hhO],[001]) plane and the one obtained in the kagome 
antifcrromagnet is evidence that these two lattices are 
very similar. This is not surprising as the kagome lat- 
tice is a cut of the pyrochlore lattice perpendicular to 
the (111) axis. As previously noted in Ref. |l] there are 
two magnetic modes coexisting in this system (see Fig. 
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I and Fig. |). The one related to the ([hhO],[001]) plane 
(Qi = [220]±[||0] or Q[ = [002]±[||0]) would coincide 
in an ordered phase with a structure where consecutive 
tetrahedra are in phase, while the one in the ([hOO],[010]) 
plane {Q2 — [210]) correspond to a tt phase between 
tetrahedra. Whereas in the quantum approach]^ there 
is a difference of amplitude between these two modes, 
within our approach both have the same amplitude so 
that the system is unable to select a particular one. The 
powder average of the neutron cross section (da/dfl), i.e., 
the average of |l^ over the directions of q is shown in Fig. 
^. The continuous line is obtained from the D — 00 
component classical antiferromagnet approach while the 
dots correspond to Monte Carlo simulations performed 
m Ref.|^. There is a good agreement between these two 
methods and they do not show any sign of magnetic 
transition. We have calculated the effect of the tem- 
perature on the powder average cross section. In Fig. 

we have plotted the variation of the profiles with re- 
duced temperature [9 = T/Tmfa where Tmfa = 2J/3). 
One can see that at temperature below T^^^, i.e., be- 
low 6 = 1, these profiles have approximatively the same 
shape. In order to control that quantity, we have rescaled 
each profile and plotted them one onto the others. As 
they almost coincide, it means that at low temperature 
we can decouple the q dependance and the 9 depen- 
dance so that {da/dn){Q,9) « {da /dVl){Q,9 = Q)f{9)_ 
The same type of decoupling was assumed in Ref. |2l| 
where the authors proposed that the spatial and the time 
dependance could be independent at low temperatures 
1^) ~ f{q)g{uj)). These properties are original and 
certainly come from the high frustration of the pyrochlore 
lattice and its related dynamics. At first approximation, 
we can state that spatial and time as well as spatial and 
thermal dependance are decoupled at energy scales lower 
than the mean field temperature. 

IV. DISCUSSION 

In this article, we have presented the structure of the 
neutron diffraction cross sections obtained from the ex- 
act solution for the D = 00 component classical antifer- 
romagnet on the pyrochlore lattice. 

In contrast to conventionnal three dimensional mag- 
nets, there is no ordering and even no extended short 
range order at T = 0. This anomalous behavior is due to 
the existence of a would be Goldstone mode at low tem- 
perature but also to the presence of two dispersionless 
excitations modes that inhibit any magnetic transition. 

The magnetic cross sections are very close to what has 
been obtained in theoretical approaches (Monte Carlo 
simulations with Heisenberg spins |39| , density matrix 
expansion for spin S = 1/2 pl[ | ) a nd experimental results 
(Y(Sc)Mn2 §, CsNiCrFe We have reproduced 

the existence of two dominant magnetic modes that could 
explain a first order transition induced by magnetoelas- 
tic effects It is remarkable that the nature of the 




FIG. 4: Neutron scattering cross section from the large-D 
pyrochlore antiferromagnet at T = in the ( [hhO] , [001] ) plane, 
(of. Fig. Dark regions are low intensity ones while white 
regions correspond to the maximum of the scattering. 

spins (classical or quantum) as well as the type of mag- 
netism (localized or itinerant) does not seem to play a 
crucial role in the pyrochlore system. Even the geom- 
etry, three dimensional here, seem to have a very little 
effect on the struture of the magnetic response. In fact, 
we obtain results similar to the ones previously observed 
on the kagome lattice. Even if this point is not solved, 
it is clear that the close relation between these two lat- 
tices, and the short range scale of the geometry are the 
relevant ingredients in the pyrochlore and the kagome 
antiferromagnets. 

Taking into account assumptions of Ref .|2^ we can say 
that spatial and thermal as well as spatial and time de- 
pendences of the magnetic cross sections are decoupled 
at low energy scale (i.e. T < J) as if S{q, T) ^ f{q)h{T) 
and S{q,uj) « f{q)g{Lu). This striking result confirms 
that this model describes an unconventionnal antiferro- 
magnet and has some similarities with a collective para- 
magnet. 

All these points strongly support that isotropic anti- 
ferromagnetic Heisenberg models (quantum or classical) 
on the pyrochlore lattice are spin liquids. 

Although this point is now clearer, it is still an open 
problem to understand the experimental results that sug- 
gest the appearance of a topological spin glass like behav- 
ior. There must certainly be other interactions or pro- 
cesses involved in this case as up to now, no sign of glassi- 
ness has ever been seen within theoretical approaches 
concerning the Heisenberg model on the pyrochlore lat- 
tice. 
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FIG. 5: Neutron scattering cross section from the large-D 
pyrochlore antiferromagnet at T = in the ( [hOO] , [010] ) plane, 
(cf. Fig. ^). Dark regions are low intensity ones while white 
regions correspond to the maximum of the scattering. 
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